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Kim and Chan have recently observed Non-Classical Rotational Inertia (NCRI) for solid 4 He in 
Vycor glass, gold film, and bulk. Their low T value of the superfluid fraction, p s /p ~ 0.015, is con- 
sistent with what is known of the atomic derealization in this quantum solid. By including a lattice 
mass density pl distinct from the normal fluid density p n , we argue that p s (T) ~ Ps(0) — p n (T), and 
we develop a model for the normal fluid density p n with contributions from longitudinal phonons 
and "defectons" (which dominate). The Bose- Einstein Condensation (BEC) and macroscopic phase 
inferred from NCRI implies quantum vortex lines and quantum vortex rings, which may explain the 
unusually low critical velocity and certain hysteretic phenomena. 
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Introduction. In 1969, Andreev and Lifshitz pro- 
posed that solid 4 He could display unusual lattice flow 
properties if the crystal vacancies and defects become mo- 
bile enough for quantum diffusion to become more signif- 
icant than thermal diffusion X They suggested that exter- 
nal objects placed within the solid could become unusu- 
ally mobile, and they developed a two-fluid hydrodynam- 
ics for the system. In 1970, motivated by Monte Carlo 
studies of the wavefunction for solid 4 He, Chester sug- 
gested that solid 4 He might exist in a periodic state that 
is subject to Bose-Einstein Condensation (BEC) and, fur- 
ther, that an alternative way to obtain BEC would be 
for the system to have vacancies in the ground stated 
Leggett shortly proposed that solid 4 He might display 
Non-Classical Rotational Inertia (NCRI), where under 
rotation the angular momentum would be less than for 
a classical stated In the rotating frame the lattice is at 
rest, but there is mass flow due to gradients in the phase 
4> of the wavefunction, which one can call "phase flow" 
to distinguish it from lattice flow. Ref. y developed a 
theoretical expression for the excess energy caused by ro- 
tation. It also indicated how to obtain an upper limit 
- and for BEC perhaps the exact value - for the super- 
fluid fraction p s /p at temperature T = 0, where p s is 
the superfluid mass density and p is the total mass den- 
sity. Estimates based on analogies to tunnelling in solid 
3 He led to values of the superfluid fraction from 0.01% 3 
to 0.0001%) 4 « but actual calculations based on Leggett's 
theory led to values from 5% to 30%A& 

Since these suggestions there have been many searches 
for "supersolid" behavior^ To date, no experiment has 
found evidence for lattice flow>2*?« However, using the tor- 
sion oscillator technique, Kim and Chan recently have ob- 
served NCRI for solid 4 He in Vycor glass. 10 Even more 
recently they have found NCRI for solid 4 He in porous 
Au and in bulkii They find that the zero temperature 
value for p s /p is on the order of 1-2%. They also find 
that the transition from the supersolid to the ordinary 
solid (which does not display NCRI) is around 0.2 K. 

The present work: (1) obtains converged values for 
the superfluid fraction at T = in terms of the localiza- 
tion of the density about the lattice sites, and using the 



experimental superfluid fraction finds reasonable agree- 
ment with independent measures of that localization; (2) 
presents a theory for the temperature-dependence of p s in 
terms of a lattice mass density pl and the temperature- 
dependent normal fluid density p n , which qualitatively 
explains how the transition temperature can be so much 
lower than in liquid 3 He; (3) proposes simple forms for 
the phonon and "defecton" contributions to p n , which 
leads to semi-quantitative agreement with experiment 
and with the observed transition temperature T c ; (4) 
proposes that quantum vortex lines can occur, and that 
they have low enough energies to explain the factor of 
700 lower critical velocities observed in the supersolid as 
opposed to the superfluidji^ (5) proposes that quantum 
vortex rings can occur, and that they may explain the 
hysteresis observed on cooling but not on heating; (6) 
proposes that adding 3 He not only lowers the superfluid 
fraction by increasing the normal fluid component, but 
also may raise T c by providing pinning sites that prevent 
vortex binding and unbinding, and thus inhibit collec- 
tive processes responsible for the transition; (7) suggests 
that quantum vortex rings in the supersolid might be 
produced, as in the superfluid, with ion sources. 

Superfluid Density at T = 0. Consider the flow 
energy E = (1/2) J p{r)v 2 s {r)dr, where p(r) and v s (r) = 

{h/m,i)V<j){r) are the local values of the mass density 
and superfluid velocity. 3 (Here is the bare 4 He mass.) 
Minimization with respect to the phase (f> leads to the 
superflow condition V • (p(r)v s (fj) = From a knowl- 
edge of p(r) and the spatial average superfluid velocity 
v s , one can determine the full superfluid velocity profile 
v s (r), and from that the flow energy E. The superfluid 
density then follows on setting E/V — (l/2)p s Vg, where 
V is the system volume. Physically, the minimization 
condition causes phase flow (from v s (r)) to occur prefer- 
entially in regions of low density, and the more localized 
the p(r), the more suppressed is p s /p from unity. Note 
that only p(f) is defined on a microscopic level; its spa- 
tially average value, the total mass density p, as well as 
p Sl are defined only on the macroscopic level. 

Ref. computed p s /p as a function of b/a, for a p(f) 
given as a sum over Gaussians of width b placed on the 
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FIG. 1: Superfluid fraction p s /p vs b/a. 

sites of a fictitious face-centered-cubic lattice of solid 4 He, 
with lattice constant a. In Refs. IT(iknd [Til the system 
was likely a hexagonal close-packed crystal or polycrystal. 
Since the hep and fee lattices have the same local packing, 
this use of an fee lattice is appropriated 

At the lower values of b/a the results of Ref. were only 
an upper bound, because the computer then available did 
not have the capability to solve for the superfiow pattern 
accurately^ Current computers permit convergent calcu- 
lations. Figure 1 gives values for p s /p for the full range 
of b/a. It gives the expected result that p s /p approaches 
unity for a delocalized system (large b/a); for b/a = 0.30, 
Ps/ p « 0.9860. It also gives the expected result that p s /p 
approaches zero for a localized system (small b/a); for 
b/a = 0.10, p s /p » 0.0009. For 0.10 < b/a < 0.20, p s /p 
is very sensitive to b/a, as shown by the inset. Over 6000 
plane waves were needed to obtain p s /p at the smaller 
values of b/a. Although p s can be a tensor for a general 
lattice, it is diagonal for a cubic lattice; for the fee lattice 
considered her, on convergence p s was independent of the 
direction of superfiow. 

The exchange considerations of Ref. give b/a near 
0.125. The inset yields that p s /p = 0.022 for b/a = 0.125. 
Considering the sensitivity of p s / ' p to b/a, this is in good 
agreement with the experimental value p s /p ~ 0.015. 
Because p s / p is so sensitive to b/a, it would be desireable 
to have a more accurate p(f) from calculations on solid 
4 He using wavefunctions possessing BECiA&i^ 

If the p(r) for finite T is employed in the above pro- 
cedure, one obtains a quantity p*(T), with p*(T — 0) = 
p s (T = 0), whose significance we discuss next. 

Superfluid Density at Finite T; Lattice Density. 
Following Ref. [lj, p s and p n are (implicitly) tensors. In 
Refll5l the author attempted to extend the hydrodynam- 
ics of RefQ to include the fact that superfiow is defined 
relative to the lattice. That work, however, did not go 
far enough because it did not explicitly include the lat- 
tice velocity vl or the lattice mass density pl in the mo- 
mentum density. To incorporate this dependence of the 



superfiow on motion relative to the lattice, we take the 
total momentum density to have the form 

g = p s v s +p n v n + p L v L = Ps(v s -v L ) + Pn(v n -v L )+pv L , 

(1) 

where Vl = du/dt is the lattice velocity and u is the 
lattice displacement relative to a fiducial state. Here the 
(net) average mass density satisfies 

p = p s + p n + p L , (2) 

and g — ► g + pv under a Galilean boost of all three 
velocities by v. Thus we take the system to initially 
possess three independent mass densities and velocities, 
rather than the usual two independent mass densities 
and velocities. Although the macroscopic dynamical 
equation a 1 ! 6 ! 16 predict that v n w vl, to develop a the- 
ory for p s and p n it is essential at the outset not to make 
this simplification. On applying the macroscopic dynam- 
ical equations, the normal modes become essentially the 
same as in RefQ, subject to the comments of LiuiiSiii 

With p* s (T) of the previous section, we identify the 
lattice density p L (T) = p(T) - p*(T). Then, by ©, 

p s (T) = p(T) - p L (T) - p n (T) = p\ (T) - Pn (T) , (3) 

rather than p s (T) = p(T) — p n (T), as holds for bulk su- 
perfluid. Hence, if p s (T c ) = defines T c , and p s /p = 
0.015 at T = 0, in the supersolid one needs only a nor- 
mal fluid fraction of 0.015 (rather than unity in the su- 
perfluid) to suppress superfiow. This would explain why 
T c for bulk solid 4 He is only about 0.2 K, whereas it is 
2.176 K for bulk liquid 4 He. In practice, at temperatures 
low enough for NCRI to be observed, p(r) may be es- 
sentially independent of temperature, so that we may be 
able to employ p*(T) w p'*(0) = p s (0) in (3). Like p and 
p s , p n arid pl are defined only on the macroscopic level. 

By contrast, consider liquid 4 He in a superleak. For 
pores wide enough that the transition is not suppressed, 
the pore superfluid density p s is geometrically sup- 
pressed but otherwise takes on bulk values, so pf^ — 
gVps, where g < 1 is a geometrical constant reflecting the 
tortuous flow pattern in the superleak, V is the porosity, 
and p s = p — p n is the temperature-dependent bulk su- 
perfluid density (so that the transition temperature takes 
on the bulk value) i& On the other hand, for pores nar- 
row enough to suppress T c , pf should reflect both the 
geometrical constraints of the geometry via g and V , as 
well as a factor that differs from p s , to reflect the micro- 
scopic geometrical suppression of the overall superfluid 
density (e.g., due to suppression or modification of the 
order parameter) m& Neither of these types of superleak 
has the same behavior as proposed in © for solid 4 Hc. 

Normal Fluid Density. As a first pass at a concrete 
expression for p n , we ignore the distinction between mo- 
mentum and crystal momentum; a proper theory thus 
must be more complex than in Ref. [HjiSi We expect 
two contributions to p n : first, from phonons (here we 
consider only longitudinal phonons, with characteristic 
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velocity w 500 m/s); and second from defect exci- 
tations, or "defectons" (perhaps "vacancions" , perhaps 
"interstitial-ons" , perhaps excitons consisting of coupled 
motions of both). We take their excitation spectrum to 
have the form 1 ; 21 ! 22 ! 23 



where A4 is the minimum energy of excitation, p is crys- 
tal momentum, and /14 is their effective mass. From 
Landaui^iSi these two types of boson excitations give 

The phonon contribution to p n is 10 -8 of p for T = 
0.2 K, and thus is negligible; the exponential term due to 
"defectons" dominates. A 4 = 2 K and m of about 60% 
of m4 give a qualitative fit to the data, thus support- 
ing this view of the origin of the normal fluid, as well as 
what is currently known about defect excitations. Tak- 
ing A4 — 1 K 25 requires a much lower effective mass to 
give a similar quality fit the data. Although JSJ explains 
the sudden onset of p n , it does not explain the long tail 
observed at higher temperatures>i2iii 

The addition to solid 4 He of a small concentration c of 
3 He should change p and p^ in proportion to c, or 

p s {T = } c) = p s {T = 0) - ape, (6) 

where a is a dimensionless constant. Moreover, at low c, 
p n likely will be suppressed in proportion to c, because 
the system has less 4 He to excite. The data have this 
qualitative behaviorplSiii 

Vortex Lines and Vortex Rings. One striking 
difference between superflow in liquid and solid 4 He is 
that the critical velocity v c is much lower for the solid 
(« 300/im/s, as opposed to rs 200 nmi/s) pi In torsion 
oscillator experiments, during an oscillation the velocity 
amplitude ranges from zero to a maximum value v m . If 
u m < v c , there is no degradation of the superflow, but if 
u m > v c , there is degradation during that part of the os- 
cillation where v > v c . The observed v c is on the order of 
what one would obtain for a single unit of circulation. 11 

The observation of NCRI implies BEC, making the 
phase cj> of the wavefunction a well-defined quantity and 
permitting the superflow pattern calculations described 
above. Another aspect of having phase as a well-defined 
quantity is that we can expect other excitations involving 
the phase. In particular, there should be quantum vortex 
lines and quantum vortex ringSp&Si This is relevant to 
the low observed values of the critical velocity. 

Just as the non-uniform density of the solid makes the 
superflow pattern non-uniform, thereby decreasing the 
average superfluid density (as in Figure 1), so the non- 
uniform density of the solid makes it easier for vortices 
to enter the system and destroy the superflow. This is 
because there are low-density regions in the solid (i.e. 



between the crystal sites), at which quantum vortices 
will preferentially nucleate .SS* 2 *! This preference occurs 
because the largest energy cost for a quantum vortex in 
a finite system is associated with the formation of its 
core, where the local density tends to go to zero. If the 
local density is already suppressed by 90% between the 
crystal sites, it costs the solid much less energy than the 
liquid to further reduce the local density to zero. With- 
out calculating in detail the vortex flow velocity field (by 
analogy to that of Ref|(J for bulk flow), it is not possi- 
ble to say how much less energy would be required to 
create a vortex line in supersolid 4 He relative to super- 
fluid 4 Hep22, However, since the energy to create a vortex 
line is proportional to p s , which at low temperature is 
suppressed from p by almost a factor of one hundred, 
and since the critical velocity is proportional to the en- 
ergy to create a vortex line^i the suppression of p s by 
non-uniform superflow provides a mechanism for under- 
standing most of the factor of 700 difference between the 
critical velocities in the supersolid and the superfluid. 11 
Presently, spread in the data due to nonuniformity of the 
He crystals prevents a meaningful comparison between 
p s and z; c i 31 i 32 i 33 

Indirect evidence for vorticity (lines and rings) comes 
from the asymmetry under changes of temperature of the 
equilibration time r between torsion oscillator runs. On 
warming to T, the r to the next run is the same as be- 
tween consecutive runs at T; but on cooling, r is larger 
than between consecutive runs at T; moreover, this time 
lengthens with decreasing temperatureji* This is consis- 
tent with thermal vorticity being produced on warming, 
but on cooling excess thermal vorticity requiring addi- 
tional time to equilibrate at the lower temperature. Some 
of this vorticity may be induced by the oscillation itself, 
in which case small amplitude oscillations will equilibrate 
faster on cooling. (Some equilibration of the system is 
due simply to the large Q ss 10 6 of the oscillator, which 
coupled with a period P m 10~ 3 s, gives a 1000 s ring- 
down time pi*) 

Although quantum vortices are well known to occur 
in superconductors^ which are solids, the vortices oc- 
cur in the phase 4> of the wavefunction for the electron 
fluid within a solid structure determined by atoms. The 
present case is more extreme, for the vortices occur in 
the 4 He, which is also responsible for the solid structure: 
supersolid 4 He may be thought of as a localized fluid, 
rather than merely as a solid. 

Effect of 3 He on the Transition Temperature. If 
quantum vortex lines and quantum vortex rings are re- 
sponsible for the details of the transition from supersolid 
to ordinary solid at about 0.2 K, then the addition of 3 He 
will provide pinning centers that suppress vortex motion 
and prevent vortex proliferation. This should make it 
more difficult to destroy superfluidity, and thus raise the 
transition temperature, as observedplSiii Such enhance- 
ment is in contrast to what happens at low tempera- 
tures, where the 3 He suppresses the superfluid fraction. 
The data also shows a long, drawn-out persistence of su- 
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perfluidity at higher temperatures ^iMi this might be a 
dynamical effect^ disappearing for an oscillator with a 
significantly longer period, for which the vortices would 
have more time to escape the pinning centers^ 

Final Remarks. Atoms in optical lattices also are 
candidates for supersolid behaviori2£'22i22iiii Their low 
moment of inertia and the imposed symmetry-breaking 
optical lattice may make it more appropriate to study 
their normal modes rather than their non-classical mo- 
ment of inertia. Nevertheless, the transition from the 
(p s = 0) Mott insulating state to the (p s / 0) BEC 
state, which can in principle be produced by decreas- 
ing the strength of the laser field producing the lattice, 
should provide a dramatic change in the moment of iner- 
tia. Likewise, in the BEC state, increasing the laser field 
strength should further localize the system, and decrease 
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Although the unusual solid flow (as opposed to "phase 
flow" ) effects suggested by Ref. 0] have not been observed, 
with p n — > p n + pl— this work may correctly provide the 
hydrodynamics of supersolids. 

It would be of great interest to find more direct evi- 
dence for the vortex excitations of the macroscopic phase 
<f>, whose existence is inferred from the observation of 
NCRI. 10 ' 11 Perhaps vortex rings can be observed using a 
method similar to that of Rayfield and Reif4i 
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